Shapefinders: a new shape diagnostic for large scale structure. 



Varun Sahni a , B.S. Sathyaprakash b , Sergei F. Shandarin 
^nter-University Centre for Astronomy & Astrophysics, Post Bag 4, Pune 411007, India 
department of Physics and Astronomy, UWCC, Cardiff, CF2 3YB, U.K. 
c Department of Physics and Astronomy, University of Kansas, Lawrence, KS 66045 



00 

On 
On 

G 



> 
m 
in 
o 

o 

oo 
On 

Pl,: 

6 

CO 



X 



ABSTRACT 

We construct a set of shape-finders which determine shapes of compact surfaces (iso- 
density surfaces in galaxy surveys or N-body simulations) without fitting them to ellipsoidal 
configurations as done earlier. The new indicators arise from simple, geometrical considerations 
and are derived from fundamental properties of a surface such as its volume, surface area, 
integrated mean curvature and connectivity characterized by the Genus. These 'Shapefinders' 
could be used to diagnose the presence of filaments, pancakes and ribbons in large scale 
structure. Their lower-dimensional generalization may be useful for the study of two-dimensional 
distributions such as temperature maps of the Cosmic Microwave Background. 



Subject headings: Cosmology — galaxies: clustering- 
methods: analytical, numerical . 



large scale structure of the universe- 



The large scale structure of the Universe is remarkably rich in visual texture. At different density 
thresholds the clustering pattern has been variously described as 'meatball-like', 'sponge-like', 'bubble-like', 



network of surfaces' etc. (Zeldovich, Einasto & Shandarin 1982, Melott 1990, ie Lapparcnt, Gcllcr 



fe Huchra 1991). Attempts to quantify this pattern, in red-shift surveys of galaxies and in N-body 
simulations, have been made using topological discriminators such as the genus curve and percolation 
statistics ( Zeldovich 1982], shandarin 1983, Gott, Melott, & Dickinson 1986; for a recent discussion see 
Sahni, Sathyaprakash & Shandarin 1997), and also by applying minimal spanning trees and statistics 
sensitive to 'shape' ( Babul fc Starkman 1992 , Luo fc Vishniac 1995 , Sathyaprakash, Sahni fc Shandarin 



1996a, Dave et al. 1997, for a review see Sahni & Coles 1995) 



Recently Minkowski functionals have been used to characterize the global geometrical and topological 



properties of pointwise distributions, e.g. galaxy or cluster catalogues ( Mecke et al. 1994 ). In this paper 
we suggest to use a set of geometrical parameters derived from Minkowski functionals to describe the 
geometry and topology of individual objects such as superclusters of galaxies. Minkowski functionals have a 
very convenient property of additivity and therefore they can be used for both an isolated structure (such 
as a cluster or supercluster of galaxies) or for a group of structures (such as all structures in an N-body 
simulation or a galaxy catalog). Probably using the Minkowski functionals themselves is the most general 
and universal approach to quantitatively representing the geometry of superclusters and voids, but often we 
are interested in the characteristic dimensions of an object as well. We have therefore devised Shapefinders 
- statistics having dimensions of [Length] as well as dimensionless statistics to characterize the morphology 
of large scale structure. 

The morphology of superclusters and voids is likely to differ for different scenarios of structure formation 
and a study of supercluster-void shapes could help distinguish between radically different alternatives 
such as gravitational instability, seed models of structure formation and models based on explosions or 
'mini-bangs'. 
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It is well known that systems evolving under gravitational instability percolate at higher density 
thresholds corresponding to progressively lower values of the filling factor, which suggests that structures in 
the percolating phase are more likely to be sheet- or filament-like since sheets and filaments occupy a larger 



surface area than a sphere at a given volume, and therefore percolate more easily (Klypin & Shandarin 



1993, Sathyaprakash, Sahni & Shandarin 1996a) For a fluid that has evolved as a result of gravitational 



instability, a low filling factor at percolation is also suggested by the Zeldovich approximation which predicts 



that the first singularities to form are pancake-like (Zeldovich 1982, Shandarin & Zeldovich 1989). However 



it is unlikely that these pancakes will be strictly planar objects. Instead it is more natural that they will 



resemble, in a manner of speaking, the curved two-dimensional surface of a cup (Arnol'd, Shandarin, & 



Zeldovich 1982, Shandarin et al. 1995). Recent work suggests that soon after pancake formation, the 



density distribution becomes dominated by filaments which act as bridges connecting neighboring clusters 



(Bond, Kofman & Pogosyan 1996, Sathyaprakash, Sahni & Shandarin 1996a), with pancakes remaining 
statistically significant^. 

Most galaxy catalogues reveal structures with typical scales ~ 50 Mpc, some such as the Great Wall, 
appear to be even bigger. The present finite size of surveys, together with the fact that most of them are 
limited to surveying galaxies within a wedge shaped region, prevents us from establishing whether the visual 
structures we see are truly filamentary (one-dimensional) or they appear filamentary because the geometry 
of the survey prevents us from acquiring a fully three-dimensional perspective (filaments in a wedge type 
survey could, for instance, be slices of two-dimensional 'sheets'). Upcoming large red-shift surveys such as 
the 2dF survey at the Anglo- Australian Telescope and the Sloan Digital Sky Survey promise to reveal large 
scale structures in their full glory and shed more light on their three-dimensional shapes. 

The importance of trying to quantify shapes of clusters and superclusters, in galaxy surveys and in 
simulations has, in recent years, led to a discussion of different statistical tools which may be sensitive to 
'shape'. While such statistics have had a measured amount of success, it is fair to say that none of them 
is entirely satisfactory. A central feature of some shape indicators is that they describe the shape of a 
collection of points (equivalently - an overdense region) by evaluating its moment of inertia tensor, which 
is similar to fitting by an ellipsoid. The ratios of the principal axes then provide a means of ascertaining 
whether the structure is oblate or prolate. This method has been widely used in determining the luminosity 
profiles of galaxies and remains a powerful tool for classifying the projected shapes of ellipticals. Its efficacy 
as a discriminator for large scale structure is, however, not quite as obvious. Whereas general physical 
principles suggest that the shapes of galaxies should be predominantly elliptical or spiral, the shape of 
overdense regions in large scale structure (clusters and superclusters of galaxies) is likely to be far more 
complicated and less prone to a classification in terms of 'eikonal' shapes such as ellipsoids. Indeed, results 
of N-body simulations show that, when viewed at different density thresholds, shapes of compact surfaces 
can vary widely, ranging from approximately ellipsoidal (at high densities), to topologically complicated 
'spongy' shapes at moderate density thresholds. 

An example of a multiply-connected surface often seen at moderate thresholds in simulations is a 



1 It is worth stressing the difference between formation of pancakes, filaments, or clumps and topological properties of density 
fields. According to the definition of Zeldovich, pancakes form in regions of three stream-flows or between shock fronts. Other 
authors however label any flattened region as a pancake. In any generic density field regions having densities above a sufficiently 
high thresholds look like compact regions (all three dimensions are similar). At a lower threshold pi the connected structure 
forms in the higher density regions p > pi, the lower density regions p < pi remain connected (in other words both percolate). 
Then, at even lower threshold p2 < pi the lower density re gions p < py cease percolating. T he second percolation transition is 



sometimes associated with the formation of pancakes (e.g. Bond, Kofman & Pogosyan 1996) 
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torus. Clearly a statistic which attempts to describe the shape of a torus by fitting with an ellipsoid would 
be widely off the mark since it would lead us to conclude that the torus has a pronounced oblate shape 
and would miss completely its tubular form - which is more like a one-dimensional filament. A tendency 
to model shapes using pre-defined 'eikonaP forms, can lead to an exaggerated emphasis of oblateness or 
sphericity over filamentarity and could easily bias our understanding of the morphology of large scale 
structure ( Sathyaprakash, Salmi, Shandarin fc Fisher 1997] ). In this paper we introduce new 'shape finders' 



which are free from the above drawbacks, and probe the shape of an object without any preordained 
reference to an eikonal shape. 

A compact surface, which could be an isodensity contour above or below a given density threshold 
in an N-body simulation, can be characterized by the following four quantities 2 : (i) Volume V, (ii) 
surface area S, (iii) integrated mean curvature: C — \ J{ki + Ki)dS, (iv) integrated Gaussian curvature 
(genus): Q = —4- J KiK 2 dS, where Ki = 1/R\ and k 2 = I/-R2 are the principal curvatures of the surface. 
Multiply-connected surfaces have Q > while simply connected have Q < 0. 

We introduce three Shapefinders each having dimensions of [Length]: Hi,i = 1, 2, 3 where Hi = V/S, 
H2 = S/C and Hz = C (for multiply-connected surfaces C/Q may be more appropriate than C). Taken 
together, the trio Hi provide robust and convenient measures of 'shape' as will be shown below. Based on 
H we can also define a pair of very useful dimensionless Shapefinders: K, = QCi, IC2) where 

Hi— Hi „ H3 — H2 
Al2 + Ati A13 + A12 

We note that < 1 by construction. JC can be regarded as a two-dimensional vector whose amplitude and 
direction determine the shape of an arbitrary three-dimensional surface. Combined with the genus, we get 
the dimensionless triad (/Ci,/C2,(5) giving information about shape as well as topology. The Shapefinders 
Hi, having dimensions of length, can be thought of as describing the spatial dimensions of an object. Thus, 
an ideal pancake (having vanishing thickness but not necessarily planar) has one characteristic dimension 
much smaller than the remaining two, so that Hi <C H2 — H3 and K, ~ (1,0). An ideal filament (a 
one-dimensional object but not necessarily straight) has two characteristic dimensions much smaller than 
the third so that Hi ~ H 2 <C H3 and K. ~ (0, 1). All three dimensions of a sphere are equal resulting 
in Hi — H2 — H3 and K, ~ (0,0). In addition, an interesting surface to consider is a 'ribbon', for which 
Hi <H 2 <H 3 and/C~ (1,1). 

The genus Q, integrated mean curvature C , and surface area S, can be derived for an arbitrary surface 
by considering its first and second fundamental forms. Consider a local coordinate patch on a surface of 
class > C 2 , r = r((9, 0), then the first (I) and second (II) fundamental forms of the surface are: 

I = dr • dr = E d9 2 + 2F dO d(j> + G d<j) 2 
II = dr • dn = L dO 2 + 2M d9d<j> + Nd<t> 2 (2) 

where n is the unit normal to the surface n = x rg / 1 x rg\, and 

E = r e ■ v g , G = • 1-0, F = i> • r e 
L = vgg ■ n, N = r>0 n, M = v^g ■ n, (3) 



= dr/d<f>, = d 2 r/d(f) 2 , etc. As a result one gets (see e.g. Lipschutz 1969 ) 

S = J J VEG - F 2 d<pd6 



2 These are also known as Minkowski functionals and have been applied to Cosmology by 



Mecke et al. 1994 . 



-4- 



KlK 2 



EN + GL- 2FM 
EG — F 2 ' 
LN - M 2 



EG-F 2 



C = J J ^±^dS, 

9 = / J Kin 2 dS. (4) 

We elaborate on the shape statistics by applying them to two surfaces - an ellipsoid and a torus. Consider 
first the triaxial ellipsoid having volume V = ^j-abc and parametric form 

r = a(sin 9 cos <f)x + 6(sin 8 sin <j))y + c(cos 9)z (5) 

where < 4> < 2tt, < 9 < n. 

Results for the trio of dimensionful Shapefinders Hi = (V/S, S/C, C) and the dimensionless 
Shapefinders JC = (JC\,K.2) are shown in Tables |l| & |^ for the triaxial ellipsoid with axis (a, b, c) and its 
deformations into a pancake, filament and sphere respectively. In all cases the Shapefinders have been 
normalized to give Jii = R (JCi = 0) for a sphere of radius R. 

Our results for idealized 'eikonaP surfaces are as follows: pancakes (a, b — a, c), a ^ c, JC ~ (1,0); 
filaments (a, b, c = b), a b, K ~ (0, 1); ribbons (a, d,c), «»()>c, K~ (1, 1); spheres (a, 6 = a, c = a), 
K, ~ (0,0). Thus, the deformation of a sphere into a pancake or filament results in the accompanying 
increase in JC: (0,0) — > (1,0) (pancake), or (0,0) — > (0,1) (filament). While deforming a filament into a 
pancake (and vice versa) it is possible to encounter an intermediate surface which resembles a ribbon, for 
this surface K, ~ (a, a), a < 1. Tables ^ & || also show an interesting symmetry between oblate and prolate 
ellipsoids. Consider an oblate ellipsoid (a, a, c), with /C = (/Ci,/C2), then, for the prolate ellipsoid (a, c, c), 
K, = (IC2, /Ci)! For shapes close to spherical, both /Ci, /C2 are small, and it is more appropriate to study the 
ratio K-i/ICj to assess departure from sphericity. Thus for pancakes /C1//C2 > 1, filaments K-2/IC1 > 1 and 
for ribbons K\/K,2 — 1. In addition the statistic 1C\ 2 , P < 1, will also accentuate small departures from 
sphericity. 

Next we consider a torus with elliptical cross-section having the parametric form 

r = (b + csin</>) cos (9 x + (b + cs'm(f>) sin 9 y + a(cos(f>)z (6) 

where a, c < b, < 4>, 9 < 2n. The toroidal tube with diameter 2nb has an elliptical cross-section, a & c 
being respectively radii of curvature perpendicular and parallel to the plane of the torus. The usual circular 
torus is given by a = c. 

Table || shows 1C & TL for a torus with an elliptical cross-section which is an interesting surface since, 
like the triaxial ellipsoid, it has three sets of numbers defining its shape but it also has a hole which makes 
it topologically non-trivial. Such surfaces (with holes) generically arise in large scale structure surveys and 
in N-body simulations, especially at moderate density thresholds. Our results for the torus are qualitatively 
similar to those for the ellipsoid but with interesting new features. For instance one has two possibilities for 
a ribbon: Ribbonl resembles a finger-ring, whereas Ribbon2 is a disc-with-large-hole, see Fig.l. (In the 
extreme case c ~ b, Ribbon2 becomes like Pancake2, and, for a ~ b, Ribbonl becomes like Pancakel.) 

Finally, for both the ellipsoid and the torus, JC remains invariant under a similarity transformation: 
(a, 6, c) — > (aa, ab, ac). 
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It is quite remarkable that TL and K, give compatible results for similar shapes obtained by deformating 
two very dissimilar bodies, a torus and an ellipsoid. We feel this indicates that the Shapefinders introduced 
by us are robust and can be used for identifying topologically complicated shapes. In a companion paper 
we shall apply the shape-statistic to N-body simulations ( Shandarin, Sathyaprakash fc Salmi 1997 )^. As 



emphasised by Sathyaprakash, Sahni & Shandarin 1996a, shape diagnosis combined with percolation theory 



provide a powerful tool with which to study gravitational clustering. Percolation theory demonstrates that 
the number of clusters/ voids in a simulation always peaks at a threshold just above/below percolation. 
Thus, it makes sense to study shapes of individual clusters/ voids at this threshold since the number of 
distinct isolated surfaces (boundaries of clusters/voids) is largest. In addition, it is likely that most surfaces 
at such low thresholds will be multiply-connected with a large genus, perhaps making it appropriate to 
define the third shape statistic to be C/Q instead of C as we have been assuming so far. 

The statistic suggested by us could also be used to study more general shapes than those appearing 
in large scale structure. For instance one could use them to study the shapes of concentrated cosmic 



magnetic fields which might have important astrophysical consequences (Ryu et al. 1997). Finally, the 
two-dimensional Shapefinders Hi = S/L, H2 = L and /C = (H2 — Hi)/(H2 + Hi) (L is the circumference 
of a curve bounding a two-dimensional area S), combined with the two-dimensional genus, could prove 
useful when studying shapes and topologies of two-dimensional contours defining 'hot and cold spots' in the 
Cosmic Microwave Background, or isodensity surfaces in projection data. (Values of K, range from zero for 
a circle to unity for a filament.) 

Acknowledgments: We acknowledge stimulating discussions with Sanjeev Dhurandhar and Somak 
Raychaudhury. SFS acknowledges financial support from the NSF-EPSCoR program and NASA grant 
NAG5-4039. 



3 It is interesting that other Shapefinders besides those mentioned above, can be constructed out of the four Minkowski 
functionals V, S, C and Q. For instance the dimensionless pair A = (Ai, A2), where Ai = /V 2 — 1, A2 = (C 2 /S) 3 — 1, is also 
a good probe of pancakes, filaments and ribbons. A shares many of the features of K. (the main difference being that, unlike 
/C, A is not bounded from above Ai j2 > 0). If like /C, we normalise A so that Ai j2 = for spheres, then, for (i) pancakes 
Ai > A 2 , (ii) filaments A 2 > Ai, (iii) ribbons Ai ~ A 2 ^ 0. 
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Fig. 1. — Deformations of a torus with an elliptical cross section. 
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Table 1: Shapcfinders for a triaxial ellipsoid with axis a, b, c. (V/S, S/C, C have dimensions of length, K,\^ 
are dimensionless.) The second column (M) is a description of the morphology of the object based on its 
dimensions: P= Pancake, F= Filament, R= Ribbon, S= Sphere. 



a, b, c 


M 


(/Ci,/C 2 ) 


V/S 


S/C 


C 


(100, 100, 3 ) 


P 


(0.83, 0.10) 


5.98 


63.9 


78.6 


(100, 3, 3) 


F 


(0.10, 0.83) 


3.82 


4.70 


50.2 


(200, 20, 2) 


R 


(0.67, 0.67) 


3.94 


20.0 


101 


(100, 100, 100) 


S 


(0.99, 0.00) 


100 


100 


100 



Table 2: Deformations of a triaxial ellipsoid with axis a, b, c. 



a,b,c {K,i,K,2) 

Sphere — ► Filament 

(100, 100, 100) (0.000, 0.000) 

(100, 80, 80) (0.004, 0.005) 

(100, 50, 50) (0.028, 0.054) 

(100, 20, 20) (0.077, 0.300) 

(100, 10, 10) (0.095, 0.540) 

(100, 3, 3) (0.100, 0.830) 
Sphere — ► Pancake 

(100, 100, 100) (0.000, 0.000) 

(100, 100, 80) (0.005, 0.004) 

(100, 100, 50) (0.054, 0.028) 

(100, 100, 20) (0.300, 0.077) 

(100, 100, 10) (0.540, 0.095) 

(100, 100, 3) (0.830, 0.100) 

Pancake — > Filament 

(100, 100, 3) (0.830, 0.100) 

(100, 70, 3) (0.800, 0.130) 

(100, 30, 3) (0.650, 0.330) 

(100, 10, 3) (0.330, 0.650) 

(100, 3, 3) (0.100, 0.830) 



Table 3: Shapefinders for a torus of radius b having an elliptical cross-section with axis a, c, (b > a,c). 
(V/S, S/C, C have dimensions of length, /C1.2 are dimensionless.) The second column (M) is a description 
of the morphology of the object based on its dimensions: Pl,P2=Pancake, F=Filamcnt, Rl,R2=Ribbon, 
S=Sphere- with- hole. 



b, a, c 


M 


(/Ci, £2) 


V/S 


S/C 


C 


(100, 99, 3) 


PI 


(0.90, 0.03) 


7.05 


136 


144 


(100, 3, 99) 


P2 


(0.88, 0.20) 


7.05 


114 


173 


(100, 3, 3) 


F 


(0.14, 0.93) 


4.5 


6.00 


157 


(150, 20, 2) 


Rl 


(0.70, 0.80) 


4.64 


25.9 


235 


(150, 2, 20) 


R2 


(0.70, 0.80) 


4.64 


25.9 


235 


(20, 19, 19) 


S 


(0.14, -0.09) 


28.5 


38.0 


31.4 



